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ABSTRACT 

We present HST Faint Object Camera observations of the asteroid 1 Ceres at near-, mid-, and far-UV 
wavelengths (A — 3636, 2795, and 1621 A, respectively) obtained on 1995 June 25. The disk of Ceres 
is well-resolved for the first time, at a scale of ~ 50 km. We report the detection of a large, ~ 250 km 
diameter surface feature for which we propose the name "Piazzi" ; however it is presently uncertain if 
this feature is due to a crater, albedo variegation, or other effect. From limb fits to the images, we obtain 
semi-major and semi-minor axes of Ri — 484.8 ± 5.1 km and i?2 = 466.4 ± 5.9 km, respectively, for 
the illumination-corrected projected ellipsoid. Although albedo features are seen, they do not allow for 
a definitive determination of the rotation or pole positions of Ceres, particularly because of the sparse 
sampling (two epochs) of the 9 hour rotation period. From full-disk integrated albedo measurements, 
we find that Ceres has a red spectral slope from the mid- to near-UV, and a significant blue slope 
shortward of the mid-UV. In spite of the presence of Piazzi, we detect no significant global differences 
in the integrated albedo as a function of rotational phase for the two epochs of data we obtained. From 
Minnaert surface fits to the near- and mid-UV images, we find an unusually large Minnaert parameter 
of fc « 0.9, suggesting a more Lambertian than lunar-like surface. 

Subject headings: minor planets, asteroids — asteroids: individual (Ceres) — ultraviolet: solar system 



1. introduction 

In the latter part of the 18th century, the Titius-Bode 
law led scientists to believe there was a "missing planet" 
in the region between Mars and Jupiter. On January 1st, 
1801 — the first day of the 19th century — Giuseppe 
Piazzi serendipitously discovered the first minor planet: 
1 Ceres 

Ceres is a G-type asteroid (Tholen 1984; Barucci et al. 
1987) that orbits the Sun with a 4.4 yr period, a semi- 
major axis of 2.7 AU, and an eccentricity of 0.097. It has 
a rotational period of 9.075 hours (Lagerkvist et al. 1989). 
Lebofsky et al. (1986) reported that thermal observations 
of Ceres before and after opposition indicate that it is a 
prograde rotator (i.e., the pole is in the ecliptic North). 
Ceres has an absolute magnitude oi H = 3.32, and a mag- 
nitude slope parameter G = 0.12 ± 0.02 (Lagerkvist & 
Magnusson 1990). Ceres has measured optical colors of 
U — B — 0.43 and B — V — 0.72, and a visual geometric 
albedo - 0.10 (Tcdcsco 1989). 

^ joel@boulder.swri.edu 

^ Permanent address: Observatoire Midi-Pyrenees, 14 avenue E. Be- 
lin, F31400 Toulouse, France 



In this paper we present Hubble Space Telescope /Faint 
Object Camera {HST /FOG) observations, which are the 
first images to resolve the disk of Ceres sufficiently to de- 
tect surface features. We present the observations, an 
analysis of the photometric results, and we discuss a few 
tantalizing features suggested by the data. Though the im- 
ages are not sufficient to provide definitive new results re- 
garding such issues as pole position and composition (due 
to restricted temporal and wavelength coverage), they do 
provide a basis data set to direct further observations that 
should answer these questions. 



2. OBSERVATIONS AND DATA REDUCTION 

The observations we report here were made by HST 
with the FOC over a 5-hour period on 1995 June 25. Ta- 
ble 1 describes the observations. 

At the time of observations, the Sun-Ceres-Earth phase 
angle was a — 19.4 deg, producing an observed illuminated 
fraction of 97.2% of the disk of Ceres; the corresponding 
defect of illumination was 0.012 arcseconds. 
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Table 1 
Summary of Observations 



Data Set 


Band 


'^center 


AA 


Filters 


Start Time 


Exp. Time 


RotationaP 






(A) 


(A) 




(1995 Jun 25, UT) 


(sec) 


Phase (j) (deg) 


x2og0101t 


mid-UV 


2795 


134 


F4ND F275W F278M 


10:02:56 


716.5 


0.00 


x2og0102t 


ncar-UV 


3626 


106 


F6ND F342W FIND F372M 


10:22:47 


910.5 


14.19 


x2og0103t 


far-UV 


1621 


159 


F175W F152M 


11:29:33 


896.5 


58.26 


x2og0104t 


far-UV 


1621 


159 


F175W F152M 


11:49:28 


1292.5 


73.61 


x2og0105t 


mid-UV 


2795 


134 


F4ND F275W F278M 


13:05:58 


1016.5 


122.67 


x2og0106t 


near-UV 


3626 


106 


F6ND F342W FIND F372M 


13:30:49 


1001.5 


139.01 


x2og0107t 


far-UV 


1621 


159 


F175W F152M 


14:42:34 


896.5 


185.87 


x2og0108t 


far-UV 


1621 


159 


F175W F152M 


15:02:29 


1292.5 


201.22 



"Rotational phase is relative to the first exposure and assumes a rotational period of P = 9.075 h (Lagerkvist et al. 1989). The given values of 
the rotational phase are measured at the midpoint of each exposure. 



Ceres' heliocentric distance was r = 2.57 AU and its geo- 
centric distance was A = 2.97 AU. At this distance, one 
arcsecond equates to about 2150 km, so the 0.01435 arc- 
sec width of one FOC pixel corresponded to a pliysical 
scale of 30.9 km on the surface of Ceres. As stated in 
the HST /FOC Instrument Handbook, the FWHM of the 
FOC PSF for the filters we used in these observations is 
approximately 0.03 arcsec, giving a Raylcigh criterion res- 
olution of ^ 0.024 arcsec (the width of ^ 1.7 pixels) or 
a spatial resolution of ~ 52 km on the surface of Ceres. 
From previous estimates of its size, we expected Ceres to 
have an angular diameter of about 0.43 arcsec (the width 
of 30 FOC pixels), giving us over 700 FOC pixels covering 
the disk of Ceres. This is significantly better spatial res- 
olution and sampling than has been obtained previously 
with ground-based adaptive optics observations (Saint-Pe 
et al. 1993; Drummond et al. 1998). 

Our HST /FOG observations consisted of two sets of ex- 
posures in three wavelength bands: the far-UV (using a 
F175W+F152M filter combination; 1621 A, AA = 

159 A), the mid-UV (F275W-hF278M; Ac = 2795 A, AA = 
134 A), and near-UV (F342W-I-F372M; Ac = 3626 A, 
A A = 106 A). Neutral density filters were used in the two 
longer-wavelength observations to keep the count rates at 
reasonable levels. We selected these UV bands because 
compared to visible wavelengths they provide better reso- 
lution, ice diagnostics, and surface contrast to search for 
features. 

The FOC images, obtained with the COSTAR optical 
correction (Jedrzejewski et al. 1994), were reduced in the 
standard STScI Routine Science Data pipeline, which per- 
formed flatfielding and geometric corrections. The result- 
ing images have a field of view of 7.3 x 7.3 arcsec^ (512x512 
pixels, with each pixel 14.35 x 14.35 mas^ in size). Fig- 
ure 1 shows the four mid-UV and near-UV FOC images of 
Ceres we obtained. The far-UV images are not shown here 
because they do not contain sufficient signal for displaying 
the disk or detecting features. However, as we show later, 
there is sufhcient signal in the far-UV images to measure 
total disk-integrated counts (and thus determine the far- 
UV flux) from Ceres, which we compare to the total fluxes 
in the other bands as a measure of Ceres' UV color slope. 



3. DISCUSSION AND RESULTS 

3.1. Regarding Shape and Pole Position 

Ceres is the largest asteroid, but even after two hundred 
years of observations, its size is still in dispute. Estimates 
of the radius of Ceres have ranged from 391 km (Barnard 
1900) to 600 km (Dunham et al. 1974). Infrared Astronom- 
ical Satellite (IRAS) observations made in 1983 imply a 
radius of 457±22 km (Matson 1986). More recent observa- 
tions have considerably improved measurements to quoted 
uncertainties better than 1%, but with values differing by 
at least a few percent. The best direct size measurements 
available are the stellar occultation observations of Mil- 
lis et al. (1987) and the adaptive optics images obtained 
by Saint-Pe et al. (1993) and Drummond et al. (1998). 
Millis et al. found projected equatorial and polar radii of 
Ri = 479.6±2.4 km and i?2 = 453 .4±4.5 km, respectively, 
and an equivalent radius of i? = \/ R1R2 = 466.3 ± 2.6 km. 
Saint-Pe et al. find values that imply Ri ~ 499 km and 
i?2 ~ 469 km, giving R = 484 km. The Drummond et al. 
observations were rotationally resolved, allowing them to 
determine a fully triaxial shape with radii of a = 508 km, 
b = 473 km, and c = 445 km, with uncertainties of about 
5 km; their projected ellipse values are about 4% larger 
than those determined from the occultation observations 
of Milhs et al.. 

We measured the disk center and the projec;ted semi- 
major and semi- minor radii of Ceres by performing el- 
lipse fits to each of the four good-signal HST images; this 
method was previously used to obtain shape parameters of 
the FS'T/WFPC2 images of Vesta (Thomas et al. 1997). 
Eac;h image was scanned along columns or lines, depending 
on which is more nearly perpendicular to the limb. The 
resulting values were used to position a sharp edge along 
the scan; this was fit using the predicted disk brightness 
and a Gaussian smear from disk brightness to pixel bright- 
ness. (The predicted brightness on the disk was simply 
taken from the scan across the edge, where pixels that are 
5-7 from the initial edge are averaged for an "on disk" 
value.) The position of the sharp edge was determined to 
0.1 pixel- width by least squares matching of the predicted 
brightness along the scan to the actual one as a function of 
the modeled position of the sharp edge. For these images. 
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Fig. 1. — The mid-UV and the near-UV HST /FOC images of Ceres. Vectors show the direction of the Sun and the image orientation. The 
values for the relative rotational phase, (p, of each image were calculated assuming a synodic rotational period of P = 9.075 h (Lagerkvist et 
al. 1989). The images have been scaled such that they are all normalized to their peak pixel value (given in Table 3) for optimal contrast and 
ease of cross-compajrisons. The Piazzi feature is near the center of the first mid-UV image {<j> = 0.0 deg). 
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Table 2 

Fitted Sub-Earth and Sub-Solar Pixel Coordinates 



Data Set 


Sub-Earth [Xc,Yc]'' 


Sub-Solar [Xs.Y,]^ 


x2og0101t 
x2og0102t 
x2og0105t 
x2og0106t 


[197.80, 398. 19J 
[196.38, 396.50] 
[192.41, 390.66] 
[190.67, 389.53] 


[200.63, 393.95] 
[199.21, 392.26] 
[195.24, 386.42] 
[193.50, 385.29] 



"The Center pixel of the disk of Ceres is determined using the limb- 
fitting procedure described in the text. The coordinate system has 
an origin such that position [1,1] is the center of the lowcr-lcft pixel 
in the FOC image. The far-UV images have insufficient S/N to 
determine the disk center. 

''The sub-solar position is AX = Xs - Xc = -1-2.83, AY = -4.24 
pixel-widths from the sub-Earth coordinate. 



obtained at modest phase angle, we performed a full-disk 
scan and corrected the down-sun (terminator) position to 
a predicted limb position of the raw data, and then ellipses 
were fit. Fits of the half ellipses (illuminated limb only) 
produced similar results, but have more uncertainty due 
to the greater uncertainty of the fit center. Table 2 gives 
the sub-Earth and sub-solar pixel coordinates in the FOC 
images resulting from our fits. 

Using this method, we obtain average values for Ceres' 
illumination-corrected, projected semi-major and -minor 
radii of Ri = 484.8 ± 5.1 km and i?2 = 466.4 ± 5.9 km, re- 
spectively, and an equivalent radius of i? = 475. 5± 3.9 km. 
Since at the time of these observations the pole is within 
1 18 deg of the plane of the sky for all of the published 
pole solutions, these radii of the projected ellipsoid are 
expected to be close to the true values. Our radii values 
are consistent with, though slightly larger than (by 1 to 
2it), the sizes obtained from the occupation measurements 
by Millis et al. (1987), and are in better agreement with 
the adaptive optics values. Some of the differences could 
be due to the effect of differing rotational phase and sub- 
Earth latitude for the different observations. Observations 
taken at different epochs sec the disk cross-section at dif- 
ferent sub-Earth latitudes and longitudes, so comparisons 
are problematic since only the Drummond et al. (1998) ob- 
servations provide an estimated 3-dimensional shape for 
Ceres. However, for an object rotating this slowly, the 
equilibrium shape is expected to be a Roche ellipsoid'^ with 
axes of a = 6 > c, so the maximum projected axis is equal 
to the the equatorial axes; the minimum projected axis is 
between a and c. 

Using our projected radii we can determine the mean 
density of Ceres for a given mass. There have been several 
mass estimates for Ceres based on its perturbations on 
other asteroids (e.g., Schubart 1974; Landgraf 1988; Goffin 
1991; Sitarski & Todororovic-Juchniewicz 1995; Viateau & 
Rapaport 1998; Hilton 1999; Michalak 2000) and on Mars 
(Standish & Hellings 1989). Table 1 shown by Michalak 
(2000) provides a complete list of mass estimates of Ceres, 
showing a range of masses from 6.7 x lO^^^Al© (Schubart 

^ Note that this is in contradiction to the results of Drummond ct 
al. (1998), who find significantly differing values for all three axes. 
Their shape values would imply a lightcurvc amplitude a factor of 
two larger than that oljsorvcd; they conjecture that albedo variations 
could decrease the shape-induced lightcurve amplitude. 



1970) to 4.3 X IQ-^^Mq (Kuzmanoski 1996), though most 
of the values are within the range 4.6-5.0 x 10^^^ 
Millis et al. (1987) estimate a mean density of p = 2.7 ± 
0.1 gm cm~^ using one of the higher mass estimates of 
(5.9 ±0.3) X IQ-'^^Mq (Schubart 1974). For this mass, we 
obtain a similar density oi p = 2.6 ±0.2 gm cm~'^. Using a 
lower mass estimate for Ceres of (4.39 ±0.05) x 1Q~'^^Mq 
from Hilton (1999), our measurements imply a density p = 
1.90 ± 0.05 gm cm-3. 

The implications of these different density values are sig- 
nificant. Given the measured shape (specifically, the ob- 
served difference between the projected radii, R1 — R2), the 
higher density is consistent with the assumption that Ceres 
is homogeneous and in hydrostatic equilibrium. However, 
at the lower density value, an object with the size and 
rotation period of Ceres would have a difference in the 
semi-major and -minor axes of o — c > 40 km due to rota- 
tional flattening. The difference of the Millis et al. (1987) 
projected radii is Ri — R2 = 26.2 ± 5.1 km, and the dif- 
ference in our fit is even less: 18.4 ± 7.8 km. Therefore, 
our results appear more consistent with the higher mass 
and density values: an object with the lower mean density 
value and the size and rotation period of Ceres should have 
greater fiattening than is observed. However, if the lower 
density value is correct, then the observed lack of signif- 
icant rotational fiattening would imply that Ceres is not 
internally homogeneous. Another possible explanation for 
this (liscrepanciy between the predicted and measured flat- 
tening would be that the pole position affected our mea- 
surement of the projected radii, such that the difference 
between the true axes, a — c, is significantly larger than 
our observed projected values, i?i — R2', if that were the 
case, then Ceres could be homogeneous at the lower den- 
sity. However, as discussed below, such a pole position is 
unlikely, again making our results more consistent with a 
higher mass value. 

Because of the sparse rotational sampling of our HST 
data, we are not able to unambiguously track sufficient 
surface features to obtain a pole solution. If we assume 
that the slight flattening we find in the profile fit is due to 
rotation, then the best fit for our data put the projected 
pole angle at about 4 deg from North, with an uncertainty 
of roughly ±15 deg. We show our pole position angle and 
ellipse fit compared to those of other observations in Fig- 
ure 2. These different values for the pole position show 
a range of over 90 deg. An initial analysis by Merline et 
al. (1996) of the motion of possible surface features seen 
in the FOC images suggested that the pole position might 
be more consistent with the Johnson et al. (1983) than 
with the Saint-Pe et al. (1993), but that result could not 
be convincingly re-established in our analysis. The un- 
certainties in our ellipsoidal fit do not allow combination 
with other results to constrain the pole. For example, 
the difference in projected axes reported by of Millis et 
al. (1987) of i?i - i?2 = 26.2 km could be reduced to 
agree with our value of 18.4 km by viewing an ellipsoid 
of a = 479.6, h = 479.6, c = 453.4 km from 30 deg lati- 
tude. However, the published range of pole positions do 
not yield any viewing at greater than 18 deg latitude.* 

^ For the instance of the Saint-Pe et al. (1993) solution, Millis et 
al. (1987) would have been observing from 3 deg latitude, and this 
study would have been viewing at 10.9 deg. 




Fig. 2. — The first mid-UV image {<f> = 0.0 deg) from Figure 1 showing the shape and pole position angles determined from various sources: 
"P" = this paper, "M" = Millis et al. (1987), "J" = Johnson et al. (1983), "S" = Saint-Pe et al. (1993), and "Dn" and "Ds" are the Drummond 
et al. (1998) north and south pole positions, respectively. 
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Thus, we may only say the shape determined here, within 
our error bars, is consistent with the occultation results, 
but does not further constrain the spin pole. 



3.2. Regarding Surface Characteristics 
3.2.1. The Piazzi Feature 

Prom examination of the images in Figure 1 , no surface 
features could be easily tracked because of the large rota- 
tional phase difference between the pairs of FOC images 
within the same filter, and comparing images between dif- 
ferent filters can lead to misidentifications due to differing 
reflectances at the different wavelengths. 

However, while not trackablo. there is one noticeable 
feature we are compelled to discuss. This is a significantly 
large, roughly circular, and centrally darkened spot, which 
appears face-on in the center of the first mid-UV image (ro- 
tational phase (j> = 0.0 deg; upper left image in Figure 1), 
and covers about 40 pixels. Using these same HST /FOC 
images, this feature was independently noted by Landis et 
al. (1998). Because of the low number of counts, this is 
clearly not an artifact of nonlinear or wrapped count ef- 
fects. It also is not a blemish from one of the FOC reseau 
marks, which are 3x3 pixels in size (smaller than this 
feature); a map of the reseau marks shows that none are 
closer than ~27 pixel-widths from the center of the disk 
of the Ceres image. 

To determine the reality of this feature, we compared 
the average number of counts per pixel within the feature 
(allowing the radius and center point to vary) to the aver- 
age number of counts per pixel in an annulus outside of the 
measurement aperture. Wo also repeated the same mea- 
surements for a annular ring rather than a photometric 
aperture. The results of these measurements are shown in 
Figure 3 and give us the following parameters: (i) the cen- 
ter point has an albedo about 30% lower than in the region 
surrounding the feature, and it monotonically increases to 
the edge of the feature; (ii) the optimal diameter of the 
feature is 8 pixel-widths (consistent with our by-eye esti- 
mate of 7 pixel- widths) , or about 250 km; and (Hi) the 
center of the dark feature is roughly 0.5 and —0.5 pixel- 
widths in the X and Y directions, respectively, from the 
center of the disk image as determined from our limb fits 
described earlier. 

We also used the Kolmogorov-Smirnov (K-S) test (Press 
et al. 1992) to examine the statistical significance of this 
feature. We generated a featureless model of Ceres based 
on the Minnaert surface fit discussed in Section 3.2.3 and 
added Poisson noise to that model. We generated a cu- 
mulative distribution function based on 1000 such random 
models, which was then used to compare to our data as 
well as to 1000 more individual noise models to determine 
the distribution of K-S statistics. Comparing the FOC 
image to the models indicates a significant (^ 6(t) differ- 
ence, implying that this feature is not likely due to ran- 
dom noise. To check that this difference is not just due 
to a poor model, we repeated the same tests after mask- 
ing out the feature's pixels in both the FOC image and in 
the models; in this case the data agreed well (^ 2a) with 
the models, indicating that the model is a good match to 
the image, and only the pixels in the area of this feature 
are statistically deviant. To test that this change was not 




Radius (pixel-widths) 



Fig. 3. — The photometric profile of the "Piazzi" feature on the 
first mid-UV image (c/) = 0.0 deg) from Figure 1. The sohd fines 
show the average number of counts per pixel (thick line) and associ- 
ated photometric uncertainty (thin lines) within a circular aperture 
as a function of aperture radius. Similarly, the dashed line shows the 
average number of counts within a 1-pixel-wide photometric annu- 
lus. The horizontal dotted lines show the average number of counts 
and associated errors in a region beyond the feature. Note the an- 
nular photometry of the feature matches this background value at a 
radius of ~4 pixel-widths, which is our estimate for the mean radius 
of the feature. 



just an effect of masking out pixels from the analysis, we 
masked out other portions of the disk, but in those cases, 
the difference between the data and model remained large 
(~ 5-9cr), again indicating that only the region containing 
the feature is statistically unique. As a final check, we ran 
all the same tests on the second mid-UV image (labeled 
"0 = 122.7 deg" in Figure 1). Since that image does not 
show any obvious feature, we would expect that the model 
would fit the data in all cases (whether we were using all 
pixels, or masking out any region) better than the 5-9cr 
deviations found for the first image with the feature not 
masked out. We found that the data in the second image 
were indeed better fit to the model in all cases, typically at 
the 3-4(7 level (only one case had a 5<t deviation), but not 
as good a match as one might hope. This could indicate 
that the model does not fit the data for the second mid- 
UV image quite as well, or perhaps the feature had rotated 
to another area of the visible disk and, though not easily 
identifiable by eye, is still statistically significant. How- 
ever, these tests indicate that the models fit reasonably 
well. 

Although our analysis cannot determine the nature of 
this feature, i.e., whether it is a crater, an albedo spot, or 
something else, we believe its existence is sufficiently estab- 
lished (to be definitive, we need higher-resolution images 
with adequate time sampling to see the motion of surface 
features). With that in mind, we propose the name "Pi- 
azzi" for this feature, in honor of the discoverer of Ceres. 

3.2.2. The Geometric Albedo 

The albedo of Ceres was determined in the standard 
manner and our results calculated from the FOC data are 
given in Table 3. The geometric albedo is defined as: 



p{X,a) = 



j^obj(A) 



/a 



(1) 
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Table 3 
Summary of Photometry 



Data Set 


Band 


Exp. Time 


PHOTFLAM'^ 


Counts 


Flux'^ 




Geometric Albedo'' 






(sec) 






Peak 


Total"^ 


(erg s"^ cm"^ 






x2og0102t 


near-UV 


910.5 


9.720899x10- 


16 


79.1 


37793 


4.035x10- 


13 


0.0558 ±0.0010 


x2og0106t 


ncar-UV 


1001.5 


9.720899x10- 


15 


86.8 


41088 


3.988x10- 


13 


0.0552 ±0.0010 


x2og0101t 


inid-UV 


716.5 


3.302684x10- 


15 


28.6 


13152 


6.062x10- 


14 


0.0298 ± 0.0006 


x2og0105t 


mid-UV 


1016.5 


3.302684x10- 


15 


42.5 


18068 


5.871x10- 


14 


0.0288 ± 0.0006 


x2og0103t 


far-UV 


896.5 


6.379518x10- 


15 


5.8 


677 


4.820x10- 


15 


0.1238 ±0.0128 


x2og0104t 


far-UV 


1292.5 


6.379518x10- 


15 


4.9 


468 


2.310x10- 


15 


0.0593 ±0.0068 


x2og0107t 


far-UV 


896.5 


6.379518x10" 


15 


4.8 


575 


4.089x10- 


15 


0.1059 ±0.0091 


x2og0108t 


far-UV 


1292.5 


6.379518x10- 


15 


5.0 


560 


2.764x10- 


15 


0.0710 ± 0.0069 



" PHOTFL AM is the inverse sensitivity factor used to convert from units of (counts sec ) to units of (erg s ^ cm ^ A ^) as described in the 
HST Data Handbook. The values hsted here are the caUbrated factors obtained from the processed image headers. 
^ Total Counts arc sky-subtracted, disk-integrated totals within a 0.29 arcsec (20 pixel-width) radius aperture. 
" Flux = (PHOTFLAM) x (Total Counts) / (Exp. Time), and is the -Fobj(A) term in Equation 1. 

The albedo values are calculated using our equivalent radius oi R = 475.5 km. Quoted uncertainties represent la values were determined 
from Poisson statistics-based aperture photometry values, the uncertainty in the measured equivalent radius, and sensitivity to how the sky 
background id determined (which is particularly significaiit for the very faint far-UV images). 



where r and A are the hehoccntric (in AU) and geocentric 
(in km) distances of the object at the time of the observa- 
tion (these values for Ceres are given in Table 3), 7?obj is 
the object's radius in km (we use on our measured equiva- 
lent radius of i? = 475.5 km), _Fobj (A) is the measured flux, 
and Fq{X) is the solar flux at 1 AU over the same band- 
pass. For Fobj we use the background-subtracted disk- 
integrated flux calculated from the aperture photometry 
measured total counts listed in Table 3. The appropriate 
UV solar flux for these observations was obtained from 
the UARS satellite UV spectrometer, and the spectrum 
was modified by the FOC instrument response function 
using the SYNPHOT tasks in iraf. The parameter /„ is 
the dimensionlcss correction for solar phase angle. The 
UV phase function for Ceres is unknown; for our analysis 
here, we will use the value of the visible phase function: 
fa = 2.56, which was calculated using the equations of 
Bowell et al. (1989) with a phase angle of a = 19.4 deg 
and slope parameter of G = 0.12 (Lager kvist & Magnusson 
1990). 

The albedo uncertainties quoted in Table 3 are Itr- values 

based on the uncertainties in the equivalent radius of Ceres 
and in the aperture photometry using a 0.29 arcsec (20 pixel- 
widths) radius aperture that has been sky-subtracted us- 
ing the modal value in a sky annulus with radii of 0.43- 
0.72 arcsec (30-50 pixel-widths);^ the FOC absolute pho- 
tometric error is roughly 10-15%. 

® Because the signal in the far-UV images was insufficient for limb- 
fitting, we were not able to determine the center coordinates of the 
disk in those images. Therefore, in our calculations of the disk- 
integrated flux, we used the center coordinates from the ncar-UV 
images; the 0.29 arcsec radius photometry aperture allowed for any 
positional shift (i.e., note that the shift between contiguous mid-UV 
and near-UV images is about 2 pixels-widths, less than 0.03 arcsec, 
and the radius of Ceres is about 0.22 arcsec). Also, because of the 
very low signal-to-noise in the far-UV images, the measured flux is 
extremely sensitive to the estimated sky value, so the uncertainties 
quoted here also include the measured scatter of values using sky 
annuli of different sizes. 



The albedo values and their uncertainties are plotted 
in Figure 4. Even within the lower limit uncertainties 
shown, the albedo at the mid- and near-UV wavelengths 
display no evidence of rotational phase dependent varia- 
tion. However, the data do show that Ceres has a red 
color slope longward of the mid-UV (0.032 per 1000 A 
for A > 2800 A); Ceres has a red color slope (though 
not as steep) through the visible bands as well (Butter- 
worth & Meadows 1985). Even within the large uncer- 
tainties in the far-UV photometry, there is a significant 
blue slope (-0.051 per 1000 A for A < 2800 A) from the 
far-UV to the mid-UV, indicative of a minimum around 
2800 A. This break in slope is consistent with the differ- 
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Fig. 4. — The UV geometric albedo of Ceres determined from the 
HST observations using a phase correction factor of fa = 2.56 in 
Equation 1. The triangles (A) indicate data from the first epoch 
of observations, and the squares (□) indicate data from the second 
epoch. The number by each point indicates the rotational phase, (j), 
of that observation, la error bars are plotted representing aperture 
photometry and equivalent radius uncertainties only; for the near- 
UV and mid-UV data the error bars are comparable to the size of 
the plotted symbols. 
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ing mechanisms for formation of features as discussed by 
Gaffey et al. (1989): absorption shortward of ^ 2500 A re- 
sults from electron-hole pair formation, whereas the longer 
wavelength (~ 2500 to 3500 A) absorptions are produced 
primarily by charge-transfer mechanisms. There is no sign 
of significant absorption due to water across the disk, as 
that feature would occur in the far-UV filter, and no sig- 
nificant drop is seen in that bandpass. However, this does 
not rule out the possibility of water ice at the poles. 

Buttcrworth & Meadows (1985) analyze International 
Ultraviolet Explorer (lUE) spectra of Ceres combined with 
data at longer wavelengths to determine albedo values 
from 2100-10000 A. Their values of p w 0.05 at near-UV 
wavelengths and p ~ 0.03 at mid-UV wavelengths are in 
excellent agreement with our results. Roettgcr & Buratti 
(1994) also provide albedo measurements of Ceres from 
lUE observations in the wavelength range 2450-3150 A, 
which nicely brackets our mid-UV measurements. Their 
calculated albedo values range from p = 0.024-0.030, with 
an interpolated value oi p = 0.025 ± 0.001 at our mid- 
UV wavelength (2795 A) whereas our averaged value is 
p = 0.029 ± 0.001. The small difference could be due to 
values used for the phase angle correction (as noted above, 
we use the visible phase function correction for these UV 
data). The agreement of these lUE values and our near- 
UV and mid-UV results give us additional confidence in 
the reality of the large blue slope we find in the far-UV. 
Roettger & Buratti (1994) point out that G-type asteroids 
have higher albedos than C-type asteroids in the visible, 
but the UV albedos of the two type are similar. However, 
our HST observations also provide the first published far- 
UV data of any asteroid, so we cannot compare the far-UV 
blue slope to any other asteroid observations in order to 
understand how unique is this feature. 

3.2.3. The Minnaert Parameters 

Since the FOC images typically resolve Ceres' disk into 
over 700 sunlit pixels, we can calculate the incidence and 
emission angles for each pixel, and in principle determine 
the Minnaert parameter, k, for Ceres. In practice, the 
brightness measured at each pixel will be a combination of 
Ceres' intrinsic albedo distribution as well as the surface's 
scattering properties. We therefore compare full-disk so- 
lutions for k for each of the near- and mid-UV images, 
which should give us results that are reasonably insensi- 
tive to surface albedo variations given the large number of 
illuminated pixels and having images of Ceres at different 
rotational phases. 

The Minnaert law (Minnaert 1941) is an empirical ap- 
proximation describing scattering from a surface. It has 
the form: 

7/F = Bo/XoV-', (2) 

where jiQ and fi are cosines of the angles of emission and 
incidence, respectively, and Bq and k are the two Minnaert 
parameters (central I/F and limb darkening parameter, 
respectively) . The parameter k typically depends strongly 
on the phase angle, although for a Lambert surface k = 1 
for all phase angles. When fitting the model to the data, 
there is a correlation between k and the estimated radius 
as discussed below; we use our effective radius value of 



Table 4 

Fit Results for Minnaert Coefficients 



Data Set 


Band 


Bo 


k 


x2og0101t 


mid-UV 


22.8 ±0.9 


0.89 ±0.03 


x2og0102t 


near-UV 


65.4 ±2.5 


0.90 ±0.03 


x2og0105t 


mid-UV 


32.1 ±1.3 


0.92 ±0.04 


x2og0106t 


near-UV 


71.4 ±2.7 


0.92 ±0.03 



Note. — The correlation coefficient between Bq and k was 0.7. 
Assumed radius is R = 475.5 km = 0.22 arcsec. 



R = 475.5 km (equivalent to 15.4 pixel- widths), although 
the fact that Ceres' disk is not circular contributes to the 
uncertainty in the result. 

We performed fits to the data using two independent 
methods: the minimization used by Veverka et al. (1989) 
in their analysis of Voyager 2 observations of Uranian satel- 
lites, and fitting a linearized logarithmic version of Eq. 2.^ 
Both methods agree within the estimated 2a uncertain- 
ties; Table 4 shows the results of our fits. If we vary the 
radius by ±10 km, the fitted values of k typically vary by 
±0.04, which are on the order of the uncertainties of the 
fits listed in Table 4. The values for Bq do not appear to 
be correlated with radius. 

The limb darkening parameter, k, is typically a strong 
function of the phase angle at the instant the observa- 
tions were made, e.g., see Fig. 4 of Veverka et al. (1989). 
Ceres' phase angle was a = 19.4 deg in these observa- 
tions. Regardless of the phase angle, however, Ceres' limb 
parameter is quite high (around k = 0.9) relative to k 
values for Uranian, Galilean, and Sattirnian satellites (Bu- 
ratti 1984; Veverka et al. 1989). The simplest interpreta- 
tion is that Ceres scatters more like a Lambertian surface 
{k = 1) than a lunar-like surface {k = 0.5), which im- 
plies that much of the light we see from the surface of 
Ceres undergoes multiple reflections. Interestingly, Saint- 
Pe et al. (1993) find just the opposite in their infrared 
images: they find that the center-to-limb brightness pro- 
file in the H and K bands is better fit by a (lunar-like) 
flat disk than a Lambertian sphere. This disagreement 
could be the result of their lower resolution and/or possi- 
ble surface reflectance differences between the IR and the 
UV wavelengths. Hestroffer & Mignard (1997) also find a 
lower value for the Minnaert parameter, k = 0.61, using a 
technique of analyzing the signal modulation of Hipparcos 
Satellite data (Ceres is not resolved in their data). Al- 
though they use an average phase angle of a = 18 deg, a 
value similar to ours, they use the smaller IRAS-hased 
radius, R = 456.5 km, which would accotmt for much 
of the difference (their Ak/AR « 0.092 km^^): using 
the radii values we find from our HST data, the calcu- 
lations of Hestroffer & Mignard (1997) would yield values 
of A; = 0.74 — 0.87 ± 0.05, in better agreement with our 
results. 



^ Veverka et al. (1989) warn against linearizing Eq. 2 by taking 
its logarithm, pointing out that this method does not preserve the 
proper weighting of the data points. However, one can preserve the 
proper weighting by appropriate propagation of the errors; specifi- 
cally, if a pixel has a value of N and an error of a, then its error in 
the log version of the Minnaert law is simply a/N. 
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4. CONCLUSIONS 

Our analysis of HST/FOG near-, mid-, and far-UV im- 
ages of Ceres have produced the following results: 

1. These are the first, well-resolved (~ 50 km) images 
of Ceres. 

2. We have made a detection of an apparently large, 

250 km diameter surface feature for which we 
propose the name "Piazzi" . 

3. The semi-major and semi-minor axes are 

Ri = 484.8 ± 5.1 km and R2 = 466.4 ± 5.9 km, 
respectively, for the projected ellipsoid. 

4. Disk-integrated photometry in the three bandpasses 

show that Ceres has averaged geometric albedo 
values of: p — 0.056 in the ncar-UV, p = 0.029 in 
the mid-UV, and p = 0.090 in the far-UV. Those 
values give a red spectral slope (0.032 per 1000 A) 
from the mid- to near-UV, and a significant a 
blue slope (—0.052 per 1000 A) shortward of the 
mid-UV. The relatively high albedo value in the 
far-UV implies there is no significant absorption 
due to ice across the disk. 

5. We detect no significant global differences in 
the full-disk integrated albedo as a function of 
rotational phase for the two epochs of data we 
obtained. 

6. From Minnaert surface fits to the near- and 
mid-UV images, we find an unusually large 
Minnaert parameter of fc w 0.9, suggesting a more 
Lambertian than lunar-like surface. 

It is clear that more such imaging observations are needed 
with better sampling of the rotation period to finally re- 
solve the continuing, long-standing uncertainty in Ceres' 
pole position, and to obtain more information on the in- 
triguing Piazzi feature detected in these observations. Also, 
more sensitive observations in the far-UV should be ob- 
tained to search for surface ice on Ceres and map its distri- 
bution; the far-UV bandpass is well-matched to the strong 
water ice absorption near 1650 A as seen both in labora- 
tory data (e.g. Hudson 1971) and in the rings of Saturn 
(Wagener & Caldwell 1988). Such observations could be 
obtained with the next generation of HST imaging space 
instruments, allowing us to probe ever more deeply into 
the nature of this largest and longest known asteroid. 
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